Oscillatory driven colloidal binary mixtures: axial segregation versus laning 
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Using Brownian dynamics computer simulations we show that binary mixtures of colloids driven in 
opposite directions by an oscillating external field exhibit axial segregation in sheets perpendicular to 
the drive direction. The segregation effect is stable only in a finite window of oscillation frequencies 
and driving strengths and is taken over by lane formation in direction of the driving field if the 
driving force is increased. In the crossover regime, bands tilted relative to the drive direction are 
observed. Possible experiments to verify the axial segregation are discussed. 

PACS numbers: 82.70.Dd, 61.20.Ja 



Phase transitions in driven systems are qualitatively 
different from their equilibrium counterparts and reveal 
a wealth of novel instabilities and pattern formations in- 
duced by noncquilibrium conditions 1]. In particular, it 
is very intriguing to resolve and follow noncquilibrium 
dynamics on the particle level which is possible for gran- 
ular grains 0, 0| , mesoscopic colloidal suspensions [J| or 
dusty plasmas driven by an external field [5j . 

In the colloidal context, recent investigations have fo- 
cussed on binary mixtures which are driven by a constant 
but species-dependent force. The latter can be realized 
by gravity @, 0] or by an electric field [§| . If the drive is 
strong enough, lanes of particles driven alike are formed 
in direction of the driving field [9(. 

In this paper we simulate a model for a colloidal mix- 
ture with an time-dependent oscillatory drive [2q ]. We 
show that this induces an axial segregation, i.e. an order- 
ing in sheets normal to the driving field (rather than par- 
allel as for laning). While such an axial segregation of two 
particle species is observed in shaken or vibrated granu- 
lar systems [HI EH [H , 13 1 , it has not yet been described 
for colloidal suspensions whose dynamics is dominated by 
overdamped Brownian motion. This makes the colloidal 
dynamics different from that of granular systems where 
inelastic collisions play an important role. At fixed den- 
sity, we map out the whole non-equilibrium steady-state 
diagram over a large range of frequencies and driving 
amplitudes and find a rich phase diagram topology. By 
increasing the magnitude of the driving force or decreas- 
ing the oscillation frequency, the system gets back either 
to a disordered state or to lane formation with an inter- 
mediate regime where mixing of segregation and laning 
as well as tilted lanes occur. Decreasing the magnitude 
of the driving or increasing the oscillation frequency, on 
the other hand, always yields a disordered steady state. 
We then study the relaxational dynamics into the steady 
states and show that first laning occurs and if at all ax- 
ial segregation occurs on a larger time scale. Finally we 
discuss the anisotropy of the mean-square-displacements 
in the different steady states as dynamical consequences 
of their different structure. 



In order to simulate the dynamics of a colloidal suspen- 
sion we utilise a two dimensional (2D) Brownian dynam- 
ics (BD) simulation with neglected hydrodynamic inter- 
actions between the particles. In particular, we consider 
an equimolar binary mixture of N = 2864 hard disks 
of diameter a. In Brownian dynamics, the time evolu- 
tion of the individual colloids is governed by overdamped 
Langevin equation with thermal noise [14]. The particles 
of sort A and B are subjected to an oscillatory inversely 
phased force f (t) = —i B (t) = fo sin (uit)e x , where u is 
the driving frequency and fo the driving strength. A sin- 
gle particle will follow a noise-averaged trajectory which 
is sinusoidal in time with an amplitude fo/£,oJ. 

We simulate a quadratic system of size L/a = 75 with 
periodic boundary conditions at a total area fraction 
(\> = Nira 2 / (AL 2 ) = 0.4. As a measure for the driving 
strength, we define the Peclet number Pe = To/rd as 
the ratio between the time Try = a 2 /(AD) it takes a col- 
loid to diffuse its own radius and the time Td = a/(2vd) 
it takes to drift the same distance under the action of 
a constant external force fo. Here D — fcsT/£ denotes 
the short-time diffusion constant where ksT is the ther- 
mal energy, £ the friction constant and Vd — /o/£ is the 
maximal drift velocity. 

The initial configuration of the BD simulation is an 
equilibrated fully mixed configuration as realized in the 
absence of the drive. At time t = 0, the oscillatory force 
is turned on instantaneously and the evolution of the 
system is followed by using a finite time-step algorithm 
[1 51 ] up to about 10 5 oscillation periods. Depending on 
the Peclet number Pe and the reduced frequency wtq, 
the system ends up in different steady states. 

The nonequilibrium steady-state phase diagram is 
shown in FigQ] spanning two decades in both driving 
frequency ujtd and strength Pe. Cooresponding movies 
showing the development from the initial to the steady 
state are avaibale in (29|. Let us discuss the topology 
of the phase diagram step by step: For very low driv- 
ing strengths, the external drive is not strong enough to 
generate order in the completely mixed equilibrium state 
leading to a disordered steady-state which is depicted 
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FIG. 1: a-d BD simulation snapshots for fixed uitd = 4 but 
different Peclet numbers Pe after 10 4 periods starting from a 
fully mixed configuration. Symbols in the left corner corre- 
spond to symbols used in (e). In a the coordinate frame is 
shown and the direction of the driving field is indicated by 
the broken arrow. In a-d the length of the solid bars corre- 
spond to the amplitude of a free particle driven without noise 
in the external field. For small Pe, we observe a disordered 
state (a), for intermediate Pe colloids segregate into stripes 
oriented perpendicular or tilted to the direction of the oscil- 
lating force (b-c), on the other hand, for high Pe, lanes are 
formed parallel to the direction of the oscillatory force (d). 
Parameters are ljtd = 4 and Pe = 2, 10, 20, 110 from a to d. 
e Nonequillibrium steady state phase diagram for fixed area 
fraction <j> = 0.4. The solid line describes the simple theoreti- 
cal estimate of the disordered-to-segregated phase boundary. 



in a typical snapshot in FigQ] a. A disordered state is 
rendered as crosses in the full phase diagram shown in 
FigHk such that the region of low Peclet numbers in- 
volves a disordered steady state. As visible in FigQ] a, 
the disordered state can possess an intrinsic finite corre- 



lation length as set by colliding clusters which is larger 
than the interparticle distance but these structures do 
not span the whole simulation box. For very high fre- 
quencies, on the other hand, a free particle would just 
perform oscillations with a very small amplitude due to 
the drive. This amplitude Vd/to is indicated as an inset in 
the corresponding simulation snapshots in FigQ] a-d. At 
high frequencies, an ensemble of many particles, is just 
rattling over a small distance and is thus not perturb- 
ing much the disordered equilibrium state. Conversely, 
at very low frequencies, lane formation similar to that 
observed earlier for a static or slowly oscillating drive 
9] occurs where the direction of the lanes coincides with 
the drive direction. A corresponding simulation snapshot 
is given in FigQ] d. In a region of finite frequencies and 
Peclet numbers, bounded by disordered and laned steady 
states, axial segregation shows up, where the stripes are 
oriented perpendicular to the drive direction, see the cor- 
responding snapshot in FigQ]b. Here the particles driven 
alike perform collectively oscillations as induced by the 
external drive but collide periodically with an opposing 
band of opposite particles. In FigQ] b, a colliding inter- 
face is seen in the middle of the snapshot while voids in 
the two other interface indicate that the bands had just 
been driven away from this position. We call this stripe 
formation perpendicular to the driving field axial segre- 
gation. Finally, in an intermediate region between axial 
segregation and laning, some steady states are observed 
with tilted bands, see e.g. the snapshot in FigQ]c, where 
a general tilt angle 9 relative to the drive direction is 
observed. However, as is discussed in detail below, the 
major part of the steady states are still either segregated 
or laned. 

There are two reentrant effects of the disordered phase 
in the phase diagram shown in FigQ]e: first by increasing 
the Peclet number at fixed frequency ojtd ~ 1, the steady 
state transforms from the disordered one into the axially 
segregated one and back to the disordered one. Second, 
now for increasing frequency at fixed Peclet number Pe < 
5 ~ 10, the disordered phase is taken over by the axially 
segregated phase and then appears again. 

A simple quantitative criterion for the disordered-to- 
segregated transition at high frequencies can be derived 
by the intuitive argument that the oscillating ampli- 
tude Vd/u> must exceed the mean interparticle distance 
cr (\/ 4>max/4> — 1) in order to induce a noticable segrega- 
tion effect. The analytical estimate of this phase bound- 
ary line is shown as a solid line in FigQ]e and agrees well 
with the simulation data at intermediate and high Peclet 
numbers. 

We now turn to the dynamics of how the different 
steady states are reached after turning on the driving 
field. For this purpose we stroboscopically determine 
the location of the colloids after each drive period T 
(at a time where all external forces vanish) and cal- 
culate the partial static structure factor S(k, nT) = 
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FIG. 2: Time evolution of a typical length scale parallel (a) 
and perpendicular (b) to the direction of the oscillatory force 
lx/cr and l y /cr, as well as of the orientation 9 of a typical 
structure (c). In the BD simulations the angular frequency is 
fixed uitd = 4 and the Peclet numbers are Pe = 20, 50, 120. 
For Pe — 50, both a tilted and segregated phase are shown. 
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FIG. 3: Distribution of the orientation 9 of the steady state 
structures starting from a fully mixed configuration after a 
simulation time of 10 4 T. For every Peclet number 200 in- 
dependent simulations have been performed. The angular 
frequency is fixed to uitd = 4. 



2/N\J2iHi e lkri | 2 in a range of wavenumbers k a G 
[-k ma x,k max ] where a G {x,y} and k rnax = 2-k/cf. Here 
rj with i G {1, . . . , N/2} denote the positions of all parti- 
cles of the same species. From the structure factor we get 
the typical length scale of the system via a first moment 
defined as 



2tt J_™^\k a \S(e a k,nT)dk a 
l a {nT)~ S{e a k,nT)dk a 



(1) 



To obtain the orientation of a typical structure in the 
system with respect to the direction of the oscillatory 
force we calculate the tensor 



-k„ 



S(k, nT) [(kk)I - k ® k] dk (2) 



(where I denotes the unit tensor and <g> the dyadic prod- 
uct) and define 6 as the angle between e x and the eigen- 
vector ei corresponding to the biggest eigenvalue of S. 
The angle 9 serves as an "order parameter" of the ob- 
served structures: in fact, 6 — points to a state of lanes 
while 9 = tt/2 corresponds to axial segregated structures 
and < 9 < tt/2 indicates a tilted state. 

The time evolution of the two length scales l x and l y 
as well as the angle 9 is shown in Fig [2] at fixed driv- 
ing frequency lotd = 4 for three different Peclet numbers 
and four different end states two two of which are segre- 
gated end-state and one is laned and tilted. Strikingly, 
the system first develops structures along the drive which 
is persistent over more than a decade in time. It is only 
after this induction time that the system goes into the 
final segregated state. This clearly shows that the for- 
mation of the axially segregated steady state is a highly 



collective phenomenon which cannot by accounted for by 
a local instability analysis. Laning, on the other hand, 
can be viewed as a local instability 0, 17 1. 



Let us now discuss the tilted state in some more de- 
tail: in the range 30 < Pe < 80 some steady states pos- 
sess a non-trivial winding around our "toroidal" system 
(which is dictated by the periodic boundary conditions 
in the two dimensional system), i.e. 9 — arctan(n/m) 
with natural numbers n, m and more precisely we observe 
n/m = 1/2,1/1,2/1. We remark that similar steady 
states were also found in constant driven diffusive lattics 
gases [H, [Hj and in binary mixtures which are driven by 
nonparallel external forces ■ We observe that with in- 
creasing Pe the probability to end in a tilted steady state 
with 9 = 7r/4 increases and gets nearly equal to that of 
an axially segregated (6 — tt/2) state for Pe — 60. In the 
range 50 < Pe < 60 there exists even a finite probability 
for a third steady state 9 = arctan(2). On average the 
tilted state is reached after 10 2 — 10 3 periods and was 
always stable during further 10 5 periods. We checked 
that doubling the system size yields again a tilted state 
with same angle 9. Though this clearly shows that tilt- 
ing is persistent, we cannot exclude at this stage that the 
tilting is an artifical finite-size effect stemming from the 
toroidal boundary conditions. However, one can spec- 
ulate [18| that the lifetime of the tilted states diverges 
with system size. If so, this would be an example for 
multistability in a finite region of the Pe — uitd space as 
a consequence of nonequillibrium dynamics (2fjj . 

We finally monitor the mean-square-displacement as a 
function of time t = nT in the region of multistability 
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FIG. 4: Log-log plot of the mean square displacement A xx 
(solid line), A yy (dotted line) for fixed uitd = 4. Green and 
black lines are for Pe = 50 and orientations 6 — n/2 and 
9 = 7r/4 respectively. Red lines represent A aa for Pe = 120 
and = 0. The long-dashed blue line correspond to the mean 
square displacement of a free colloid. 



and define the corresponding displacement correlations: 



N 



A al3 (t) = l/N^2{An(t) ■ e a )(ATi(t) ■ 



(3) 



where Ar^(t) = r^(t) — rj(0) and a, (3 6 {x,y}. More 
precisely we consider systems which end up in one of 
the three steady states with 9 = 0, 7r/4, ir/2. The two 
cases 7r/4, tt/2 are calculated at the same Peclet number 
Pe = 50. Results for A a p(t) are presented in FigfJ] on 
a double-logarithmic scale. We conclude that for small 
times the motion is diffusive while for long times it is ei- 
ther diffusive or bounded. A quasi-bounded motion arises 
if the mean-square-displacement is taken for motion per- 
pendicular to the structures formed. In comparing the 
axially segregated state with the tilted state, the mean 
square displacement of the two systems is indistinguish- 
able during the first 10 2 periods, but enhanced compared 
to a free colloid due to enormous number of interparticle 
collisions which increase the particle mobility. Moreover, 
A xx > Ayy for small times due to initial local laning. For 
n > 10 2 and 9 = tt/2, particles can unhindered diffuse 
in y direction and are trapped in x direction (A xx nearly 
constant), while for 9 — it/ 4 the particles diffuse equally 
in x and y direction. Comparable to a finite mass cur- 
rent along the tilted interface in driven lattics gases [l8| , 
we notice that for 9 — ir/4 the offdiagonal part of the 
MSD does not vanish A xy ^ which is a consequence of 
a shear-flow-like displacement field. This is opposed to 
9 = 0, tt/2 where A xy vanishes. 

In conclusion, we have shown that two species of os- 
cillatory driven colloids segregate into stripes perpendic- 
ular to the drive for appropriate frequencies and driv- 



ing strengths. We have performed additional simulations 
which show that axial segregation is stable in various sit- 
uations: i) for softer interparticle interactions, ii) in three 
spatial dimensions where the stripes are two-dimensional 
sheets of finite width, iii) in slit-geometry confinement 
when the driving is perpendicular or parallel to the slit, 
iv) if hydrodynamic interactions mediated by a solvent 
flow are included in the simulation and the physical vol- 
ume fraction is small. The disordered state can be reen- 
trant for increasing the driving frequency and strength. 
For increasing strength there is a transition towards lane 
formation delineated by a mixed situation of laning and 
segregation where tilted stripes do also occur. 

In principle, it is possible to verify the segregation ef- 
fect in real-space experiments on driven binary colloidal 
mixtures. Experimental realisations are superparamag- 
netic colloids driven by a gradient in a magnetic field 
PH , oppositely charged colloidal mixtures exposed to an 
alternating electric field colloids driven by gravity 
[H, [2l[ in a rotating cell or colloidal mixtures driven by 
dielectrophoretic force in a nonuniform AC electric field 
[22I In particular we think that opposite charged 

colloids in an AC electric field are a very promising real- 
ization to see axial segregation [3fj| . Since hydrodynamic 
interactions are strongly screened in this case 24j, our 
model suitably generalized to three spatial dimensions 
and Yukawa pair interactions should be appropriate to 
describe this system. 

It would be interesting to explore higher densities 
where freezing occurs [l9| and to investigate the solvent 
flow field which might strike back to the pattern forma- 
tion at intermediate physical volume fractions. Finally, 
it is challenging to describe axial segregation by a micro- 
scopic theory based on the Smoluchowski equation where 
either dynamical density functional theories [l6l [Tm] or 
mode-coupling approaches [25j could be used. We leave 
this important study for the future. 
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